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The star-forming progenitors of massive red galaxies 
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ABSTRACT 

The link between massive red galaxies in the local Universe and star-forming galaxies 
at high rcdshift is investigated with a semi-analytic model that has proven successful 
in many ways, e.g. explaining the galaxy colour-magnitude bimodality and the stellar 
mass-age relation for red-sequence galaxies. The model is used to explore the processes 
that drive star formation in different types of galaxies as a function of stellar mass and 
rcdshift. We find that most z = 2 — 4 star-forming galaxies with M* > 10 10 Af Q evolve 
into red-sequence galaxies. Also, most of the massive galaxies on the red-sequence 
today have passed through a phase of intense star formation at z > 2. Specifically, 
~ 90% of today's red galaxies with > lO n M were fed during this phase by 
cold streams including minor mergers. Gas-rich major mergers are rare and the effects 
of merger-driven starbursts are ephemeral. On the other hand, major mergers are 
important in powering the most extreme starbursts. Gas-rich mergers also explain the 
tail of intermediate-mass red galaxies that form relatively late, after the epoch of peak 
star formation. In two thirds of the currently red galaxies that had an intense star 
formation event at z < 1, this event was triggered by a merger. 

Key words: galaxies: ellipticals — galaxies: evolution — galaxies: formation - 
galaxies: high-rcdshift — galaxies: starbursts 



1 INTRODUCTION 



Sub-mm observations have uncovered a new population of 
high- 2 star- forming galaxies (SFGs), which are detected 
through the far-infrar ed emission of dust heated by UV - 
bright newborn stars (|Smail et al.lll997l ; lEales et aT1ll999l ). 
These lie along a sequence in the star formation rate (SFR)- 
stellar m ass plane that is elevated with respect to the z = 
relation (|Noeske et alJ^OOTbtfElbaz et al.ll201ll h 

In the local Universe, ultraluminous infrared galaxies 
(ULIRGs), whose SFRs are several times higher than lo- 
cal SFGs of the same mass, are directly linked to merg- 
ers (|Sanders et al.|[l986l h Hydrod ynamic simulations repro- 
duce this observational finding jMihos fc Hernquistl 1 19941 : 



iBarnes fc Hernquistl Il991l : iMihos fc Hernquistl 1 19961 ). The 

same simulations also find that the morphological properties 
of the remnant of gas-rich mergers are consistent with thos e 
of L, ellipticals (|Barnes fc Hernquistlll996l ; ICox et ai1l2006h . 
These observational and theoretical results have been used 
to support the merger scenario, i n which elliptical galaxies 
are formed by merg ers of spirals (|Toomre fc Toomrelll972l ; 



iHopkins etlj1l2006t ). In this scenario, ULIRGs correspond 
to the formation of elliptical galaxies caught in the act. 



Massive red galaxies (i.e., Es and SOs) for med the bulk 
of their stars at high z (|Thomas et al.|[2005l h where con- 
ditions were different from the local Universe. Therefore, 
it is not straightforward to generalise the lesson from lo- 
cal ULIRGs to higher redshifts. However, progress in astro- 
nomical instrumentation has made it possible to observe the 
morphologies and kinematics of the high-z counterparts of 
ULIRGs. They are galaxies that are detected in the sub- 
mm and that form stars at rates of hundreds of solar masses 
per year. We refer to high-z ULIRGs as sub-mm galaxies 
(SMGs), independently of their actual sub-mm fluxes. 

lEngel et all (|2010h have used CO interferometric data 
to conclude that most bright (Ltr > 5 x 10 12 L W ) SMGs 
are major mergers. iForster Schreiber et all (|2009l ') have ob- 
served a larger sample of less extreme objects with inte- 
gral field spectroscopy. They have found that only one third 
of their galaxies are interacting or merging. Another one 
third are rotation-dominated turbulent d iscs and the rest 
are ve locity-dispersion-dominated objects. iKartaltepe et all 
|201ll ) present further evidence that the highest-SFR galax- 
ies tend to be highly disturbed even at z ~ 2. 

These observations have a natural explanation in a 
scenario in which cold streams in massive haloes are the 
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main mode of galaxy mass accretion (|Keres et al.l 120051 ; 
iDekel fc Birnboimll2006fh The cosmolog ical simulations and 
theoretical analysis of Dekel et al.l (|2009l ) show that only one 
third of the stream mass is in clumps leading to mergers of 
mass ratios greater than 1:10. The rest is in smoother flows 
that preserve the disc's rotation. However, discs formed by 
fast accretion of cold streams are turbulent and violently un- 
stable. They break into clumps that coalesce into a central 
spheroid, hence providing an altern ative to the merger sce - 
nario for the forma tion of spheroids (|Elmegreen et al]|2008f ) . 
IDekel et al.l i|2009l ) find that most (~ 3/4) high-z SFGs are 
stream-fed and that mergers are only necessary to explain 
the brightest SMGs, the latter being the highest-SFR sub- 
class within SFGs. 

The question is which of these two types of objects are 
the progenitors of z = massive red galaxies, and conversely, 
how many z — massive red galaxies have such progeni- 
tors. We should also like to know what proportion of high- z 
SFG's ultimately turn into massive red galaxies. In this ar- 
ticle, we address these questions with a semi-analytic model 
that has proved successful in explaining local data, such as 
the galaxy colou r - magnitude distribu tion in the Sloan Dig- 
ital Sky Survey (|Cattaneo et a l. 2006) a nd the stellar mass - 
age relation for red-sequence galaxies (|Cattaneo et al.ll2008h . 
The key feature that makes the model successful is inclusion 
of multiple mechanisms (explained below) whereby star for- 
mation can be quenched in SFGs, causing them to move to 
the red sequence. The model thus provides a plausible labo- 
ratory in which the importance of various gas accretion and 
quenching mechanisms for different types of galaxies can be 
compared with one another. 

In this paper, the model is used both to evolve the prop- 
erties of high-z SFGs to low z and to trace the progenitors of 
today's massive red galaxies back in time. It is necessary to 
analyze the problem both ways because even if all high red- 
shift galaxies of one type evolved into massive red galaxies, 
that would not imply that all massive red galaxies derived 
from that type of object. We shall also see how the model 
can be used to separate stream-fed and merger-driven star 
formation. 

It is straightforward to evolve galaxies from a given z 
to 2 = 0, once a semi-analytic model is given. It is trickier 
to tell what type of object are the progenitors of the local 
population of massive red galaxies because the properties of 
the progenitors of massive red galaxies are time-dependent. 
The main progenitor of a z — elliptical may be a clumpy 
disc at z = 4, a major gas-rich merger at z = 3, and a 
massive red object at z = 2. 

We resolve this ambiguity by introducing the redshift 
Zpcak at which the star formation rate (SFR) of a galaxy's 
main progenitor has an absolute maximum. Stream-fed star 
formation is much more continuous than merger-driven star 
formation. If it dominates at z poa k, it dominates at all times. 
In this case, star formation is unambiguously stream-fed. 
The problem is more complicated when the peak SFR is 
linked to a merger event. In that case, we need to investigate 
the event's contribution to the build-up of the final stellar 
mass. 

The plan of this work is thus as follows. In Section 2, 
we recall the main assumptions of the GallCS semi-analytic 
model used for this work. Our goal is not to give a de- 
tailed description of the GallCS model (it can be found in 



lHatton et aT]|2003l and lCattaneo et alj|2006h . It is rather to 
make clear how we separate stream-fed and merger-driven 
star formation. We shall see that, in our model, stream-fed 
gas accretion and minor merger contribute to both disc star 
formation (directly) and bulge star formation (indirectly, by 
triggering disc instabilities), while major mergers only con- 
tribute the triggering of bulge star formation. Therefore, the 
relative importance of bulge star formation can be used to 
put an upper limit to the importance of major mergers in 
the star formation histories of galaxies. In Section 3, we de- 
scribe how we assign a galaxy descendent population to a 
parent one, and vice versa. In Section 4, we present the mass 
- SFR relation at z = 2, 3, 4, and we discuss which subset 
of the galaxies that populate this relation ends up on the 
red sequence at z = 0. In Section 5, we discuss the progen- 
itors of massive red galaxies, show some characteristic SFR 
histories, and describe how the properties at 2 pB ak (z pea k, 
SFRp 0a k, driving mechanism) depend on a galaxy's stellar 
mass at z = 0. Section 6 presents the conclusions of the 
article. 



2 THE GALAXY FORMATION MODEL 

GallCS (Galaxie s In Cosmological Simulations; 
lHatton et all 120031 ') is a method to simulate the for- 
mation of galaxies in a ACDM Universe. It combines 
cosmological N-body simulations of the gravitational 
clustering of the dark matter with a semi-analytic (SAM) 
approach to the physics of the baryons (gas accretion, 
galaxy mergers, star formation and feedback). 

The version of GallCS used here is the sam e as 
the 'new mode l ' intro duced in ICattaneo et al.1 i|2006l ') and 
ICattaneo et al.l (|2008l 'l and uses the same parameter values. 
We now recall its fundamental assumptions. We concentrate 
on the points that are relevant for this work and refer to the 
articles above for a more detailed description. 



2.1 Dark-matter simulation 

The cosmological N-body simulation that follows the hierar- 
chical clustering of the dark-matter component was carried 
out with a parallel tree code. It assumes a flat ACDM Uni- 
verse with a cosmological constant of Qa = 0.667 and a 
Hubble constant of H = 66.7kms~\ The ACDM power 
spectrum of initial fluctuations is normalized to a& — 0.88. 
The computational volume is a cube of size (150 Mpc) 3 with 
256 3 particles of 8.3 x 10 9 M© each and a smoothing length 
of 29.3 kpc. The simulation produced 100 output snapshots 
spaced logarithmically in expansion factor (1 + z) _1 from 
z = 35.59 to z = 0. 

Eac h snapshot was ana lysed with a friends-of-friends al- 
gorithm (|Davis et al.ll 19851 ) to identify virialized haloes con- 
taining more than 20 particles. The minimum halo mass is 
thus 1.65 x 10 n M Q . 

Merger trees are constructed by linking the haloes iden- 
tified in each snapshot with their progenitors in the previous 
snapshot, i.e., all predecessors from which the halo has in- 
herited one or more particles. 
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2.2 Disc formation at the centre of dark-matter 
haloes 

The baryons in each halo are assumed to cool efficiently, i.e. 
to stream cold, onto the central galaxy until the dark matter 
halo reaches a critical mass of 



Merit = M shock x max{l, 10 



1.3(z 



}■ 



(1) 



This assumption has a simple physical justification. Above 
a critical halo mass of M s hock ~ 10 12 Mp), the gas that 
strea ms into the halo is shock heated (jPekel fc Birnboiml 
l2006h and becomes vulnerable to black-hole heating, which 
preve nts it from cooling do wn again (see Cattaneo et all 
20091 for a review; also s ee iDekel fc Birnboiml 120081 and 



Khochfar fc Ostrikenlioosl for a discussion of the contribu- 



tion of gravitational heating). 

The e ffects of changing M cr it have been studied quanti- 
tatively in lCattaneo et all (|2006l ) . where we have shown that 
the value of this parameter is constrained within a factor of 
two. 

The term oc 10 1,3z was introduced to account for the 
more efficient penet ration of cold streams in m assive haloes 
at hig h redshift (see lDekel fc Birnbo im 2006 and lDekel et al.l 
120091 ). 

In lCattaneo et all |2006h . we set the model parameters 
to M s hock — 2 x 1O 12 M0 by fitting the colour- magnitude 
distribution in the Sloan Digital Sky Survey and to z c = 3.2 
by fitting the Lyman-break galaxy luminosity function at 
z ~ 3. The gas that falls to the centre settles into a disc. 
Star formation is activated when the gas surface density is 
S saa > 10 20 H cm 2 ~ 1 Mq pc 2 . The star formation rate 
is described by the star formation law: 

M coW 



M 



(2) 



P* £dyn 

where M co id is the mass of cold star-forming gas and idyn 
is the half-rotation time at the disc's half-mass radius (de- 
termined assuming angular momentum conservation and an 
exponential profile). The free parameters /3 = 50 and a = 0.6 
were fixed by fitting the Kennicutt law (see lGuiderdoni et al.l 
1 19981 ) and the luminosity function of Lyman-break galaxies, 
respectively. It is important to be aware that our predic- 
tions for the luminosity function of Lyman-break galaxies 
at z ~ 3 are highly sensitive to our dust extinction cal- 
culations, which contain large uncertainties, and that these 
uncertainties trickle down in the best-fit values for both a 
and z c - 

In a recent article, iKrumholz et all (|2012T ) studied how 
Eq. which is equivalent to the Kennicutt law for local 
discs, extrapolates to high-redshift discs and starburst. They 
found that Eq. ([2]) applies to all galaxy types and at all red- 
shifts with q» ~ if tdyn is the local freefall time. However, if 
tdyn is the global dynamical time, i.e. the half-rotation time, 
then Eq. ([2]) extends to high-redshift discs, but underesti- 
mates the SFR for both low- and high-redshift starbursts. 
Therefore, the local freefall time is a better estimator of the 
star formation timescale than the global dynamical time. 
However, in our model we are obliged to use the global time 
because our model does not consider the vertical structure of 
galactic discs, which is important for the local freefall time. 

As Eq. @ with a, ~ extends reasonably well to high- 
z galaxies classified as discs, one may argue that there is no 
ground for allowing a* > 0. However, if many of the high- z 



galaxies classified as starbursts are starbursting discs (very 
much like the local galaxy M82) , then there is indeed a case 
to do so. In that case, boosting the mean star-formation ef- 
ficiency of discs at high z would be a way to account for the 
fact that, at high z, a much greater fraction of the disc pop- 
ulation is in a starbursting rather than a quiescent mode. 
Notice, however, that, for a = 0.6, the term (1 + z) a boosts 
the SFRs at z ~ 3 by little more than a factor of two, and 
that this is less than the precision with which the normali- 
sation of the Kennicutt law can be determined at high z. In 
the Discussion, we shall consider how the uncertainty on z c 
and a, may affect the conclusions of the article. 

As smooth gas accretion is the only mechanism through 
which discs are allowed to acquire gas, disc star formation is 
always stream-fed star formation. This point will be of great 
importance for the interpretation of our results. 



2.3 Bulges 

In GallCS, bulges are formed by two mechanisms: mergers 
and violent disc intabilities. 

Bulge growth via disc instabilities is modelled by as- 
suming that the bulge mass increases until it is sufficiently 
large to stabilize the disc. The stability criterion is 



0.5GM disc GM(r 1/2 ) 
- < V\ 



ri/2 



ri/2 



(3) 



where Mdisc is the disc mass, ri/2 is the exponential disc's 
half-mass radius, M(ri/ 2 ) is the total mass within the disc 
half-mass radius, and r\ is a parameter that cont rols the in- 
stability threshold (see e.g. Ivan den Bosch] 1 19981 ). This cri- 
terion translates into the condition 

V0.5M disc < 7?v/0.5Mdi S c + Mbuige + M dm (ri /2 ), (4) 

where M d m(r 1/ / 2 ) is the mass of the dark matter within the 
disc's half-m ass radius. The inst ability threshold is set to 
rj = 0.7, as in lHatton et al.l (|2003h . For this value of 77, blue 
galaxies have bulges, yet they remain largely disc-dominated 
objects. The influence of the rj parameter on our results will 
be discussed in Section 6 (Discussion and Conclusion). 

Eq. @ implies that, in the absence of mergers, the 
growth of the bulge is directly linked to the growth of the 
disc. Since bulges formed by violent disc instabilities never 
dominate the total galaxy mass in our model, bulge star 
formation triggered by disc instabilities is always secondary 
with respect to disc star formation. This is true irrespective 
of the star formation law. Even if bulge star formation hap- 
pens on a much shorter timescale than disc star formation 
(see below), gas cannot be converted into stars faster than 
it is supplied to the bulge. In a context where the bulge gas 
accretion rate is more or less proportional but always lower 
than the disc gas accretion rate, disc star formation is bound 
to be dominant. 

In GallCS, galaxy mergers are mainly due to orbital de- 
cay of satellite galaxies through dynamical friction after halo 
mergers have resulted in the formation of galaxy groups and 
clusters. The fraction of the disc that is transferred to the 
bulge in a merger grows with the mass ratio of the merg- 
ing galaxies. It ranges from zero for a very minor merger 
to unity for an equal mass merger, which means that bulge 
formation is essentially linked to major mergers. The bulge 
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that forms is assumed to have a lHernquistl <|l990f) profile and 
its radius is determined based on an energy conservation 
argument. This simplified picture of the dynamics of mor- 
phological transformations is, nevertheless, consistent with 
key observational constraints such as the Faber - Jackson rela- 
tion and the Fundamental Plane of spheroids (|Hatton et alj 
l2003h . 

Readers who are familiar with the GallCS model will 
know that, in the code, gas is not moved directly from the 
disc to the bulge but it passes through a transitional compo- 
nent called the starburst. The star formation law has exactly 
the same form (Eq. for all three components, but the star- 
burst is assumed to have scale length that is equal to one 
tenth of the disc scale-length. Therefore, the starburst dy- 
namical time that enters Eq. © is ten times shorter than the 
bulge dynamical time, which implies a factor of ten increase 
in the SFR. Stars formed in the starburst are moved to the 
bulge after they have reached an age of 100 Myr and the only 
gas in the bulge is that from stellar mass loss. However, it is 
not necessary to enter this level of detail to understand any 
of the conclusions of this article. We shall therefore simply 
speak of discs and bulges, where by bulge we mean the sum 
of the bulge and the starburst. 

It follows from our discussion that while disc star for- 
mation is always stream fed, bulge star formation can be 
caused by two processes, since both gas-rich mergers and 
violent disc instabilities (such as those deriving from fast 
accretion of cold streams) contribute to bulge star forma- 
tion. This point is important for the interpretation of our 
results and will be discussed further later on. 



2.4 Feedback 

We conclude our presentation of the GallCS model with 
a brief description of how it handles stellar evolution and 
feedback. 

Stars are evolved between snapshots using substeps of 
at most 1 Myr. During each sub-step, stars release mass and 
energy into the interstellar medium. Most of the mass comes 
from the red giant and the asymptotic giant branches of 
stellar evolution, while most of the energy comes from shocks 
due to supernova explosions. The enriched material released 
in the late stages of stellar evolution is mixed with the cold 
phase, while the energy released from supernovae is used to 
reheat the cold gas and to return it to the hot phase in the 
halo (Eq. [5]l . Reheated gas is ejected from the halo if the 
potential is shallow enough. The rate of mass loss through 
supernova-driven winds M w is determined by the equation 



1-2 

-M W V CBC — eSN77SNi?SNM s t ar , 



(5) 



where _Esn = 10 erg is the energy of a supernova, ?7sn = 
0.0093 is the number of supern ovae for 1 M Pl of sta rs formed 
and «osc is the escape velocity l|Dekel fc Silk|[l98r3 ). 

In GallCS, feedback is computed separately for each 
galaxy component. We use v csc ~ 1.84w c for discs and 
v csc = 2a for bulges/starbursts. The supernova efficiency 
esN — 0-2 is similar to that commonly adop ted in SAMs 
(ISomerville fc Prim ack 19991; ICole et alj|200oh . 



2.5 Quenching of star formation 

As in Cattaneo et al. (2008), in haloes above M cr it, we do 
not just shut down gas accretion; we also shut down star 
formation. Moreover, we suppress gas accretion in galaxies 
where the bulge mass is larger than half the total stellar 
mass. As long as Mhaio < M cr it, gas accretion can start 
again if this condition is no longer verified. 

These further assumptions were introduced not based 
on a compelling physical argument, (though some physical 
justifications are possible; see below), but rather because 
we found that they improved the agreement with the ob- 
served galaxy colour - magnitude distribution, even though 
an acceptable fit could be obtained without them (Fig. 9 of 
Cattaneo et al. 2006 shows how these further assumptions 
improve the basic model in which we simply shut down gas 
accretion at Mh a i > M cr it)- 

The second assumption implies that many (> 57%) red- 
sequence galaxies are quenched following a merger event that 
has caused the bulge mass to increase above the disc mass 
(note that the bulge mass never exceeds the disc mass by vi- 
olent instability alone; see previous discussion). This is par- 
ticularly true for galaxies around M* ~ 10 1 Mq (Cattaneo 
et al. 2008, Fig. 3). 

Two possible physical explanations for this behaviour 
are i) quenching induced by quasar fee dback following black 

120051 : 
for a 



hole growth activat ed by me rging dSpringel et al. 
iHopkins etafl 120061 : also see ICattaneo et al] 12006 
review), and ii) m orphological quenching as proposed by 
iMartig et al] (|2009l '). These authors presented a picture in 
which a large bulge stabilizes disc instabilities and therefore 
inhibits the formation of spiral arms, which are the main 
sites of star formation in a spiral galaxy. 

In the first case, star formation is quenched because the 
cold gas is blown out. In the second, the gas is not blown 
out, but it is prevented from making stars. One should no- 
tice, however, that shutting down gas accretion and shut- 
ting down star formation produce very similar effects when 
plugged into our semianalytic model. 



3 BRIDGING PAST AND PRESENT 

Having described the model, we are now ready to explain 
how we analyse the results. As this work is focussed on 
the nature of high-z SFGs and the progenitors of low-z red 
(early-type) galaxies, we need a procedure to assign z — 
descendants to high-z galaxies and high-z progenitors to 
z = galaxies. One way to do it is to follow the galaxy 
flow on a stellar mass - redshift diagram. 

Fig. [T] has been constructed by separating the stellar 
mass - redshift diagram into tiles. For each tile, that is, for 
each a stellar mass and redshift interval, we take the galax- 
ies one by one and ask what is the median mass of their 
immediate descendants. Once we have this information, we 
can use a trait to connect each tile to the tile in the following 
timestep that corresponds to the median descendant mass. 
This is what we have in the left panel of Fig. [1] which shows 
the 'forward' flow of galaxies on the M* - z diagram. 

The middle panel is very similar, but it shows the 'back- 
ward' flow. Instead of looking for the median mass of the 
descendants of the galaxies on a tile, we look for the median 
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Forward Descendants Backward Progenitors (All) Backward Progenitors (Red) 




0246 0246 0246 



z z z 

Figure 1. Forward (left panel) and backward (central and right panels) evolutionary tracks of galaxies on the M» - z plane. The left 
panel and the central panel are for the entire population. The right panel is only for galaxies that are on the red sequence at z = 0. 
Forward evolutionary tracks are constructed by linking a tile with the tile corresponding to the median mass of the descendant galaxies. 
Backward evolutionary tracks are constructed by linking a tile with the tile corresponding to the median mass of the main progenitors. 




mass of their main progenitors. So we connect each tile with 
another tile at the previous timestep. 

The right panel is identical to the middle panel, except 
that it only traces back in time the progenitors of those 
galaxies that are red at z = 0. 

Let us start by comparing the left and middle panels. 
We see that the forward and the backward evolution are not 
identical. The reason for this difference is easily explained. 
It is possible that most of the galaxies in two mass bins, Mi 
and M2, at redshift z\ come from the same mass bin, M3, 
at the previous redshift step 22 > z\ . That means that some 
of the galaxies in M3 at Z2 end up in Mi at z\ and some 
in M2. Let us assume that most end up in Mi. Then only 
Mi at z\ will be linked to M3 at Z2 in the forward evolution 
but both Mi and M% at z\ will be linked to M3 at Z2 in the 
backward evolution. 

This is what happens, for instance, in the mass bin 
M, ~ 10 11 ' 5 - 10 11 ' 55 M at z ~ 3.8 - 4. The very few 
galaxies on this tile are galaxies that belonged to the mass 
bin M* ~ 10 11 ' 3 - 10 11 ' 35 M Q at the previous timestep and 
that have experienced a sudden mass increase of ~ 60% due 
to a major merger (Fig. [1] middle panel). However, most of 
the galaxies with M. ~ 10 11 ' 3 - 10 11 ' 35 M Q at z ~ 4 - 4.2 
do not experience a merger at the next timestep (Fig. [T] left 
panel). 

Let us now compare the middle panel (the progenitors 
of all galaxies) to the right panel (the progenitors of red 
galaxies). At high M*, we hardly see any significant differ- 
ence. This should not surprise us because, at z = 0, most 



galaxies with M* > 10 11 M are red. At M* < 10 10 8 M , 
we see that most red galaxies (right panel) have not grown 
in mass for the last half the cosmic lifetime (since z ~ 1). 
In contrast, the overall galaxy population (middle panel) is 
dominated by objects that are still growing or that have 
grown until very recently. The same difference is also seen 
when we compare the left panel and the middle panel. 

While we have spent three paragraphs to clarify these 
differences, the main result that emerges from a comparison 
of the right and the left panel in Fig. [T] is, in fact, their simi- 
larity, particularly at high masses, which constitute the focus 
of this article. Hence, galaxy histories can be described as 
one-parameter family of curves, with the stellar mass at any 
given redshift as the key quantity that determines on which 
evolutionary tracks a galaxy lies. This is also the reason why 
we can establish a one-to-one relation between z = galaxies 
and their high-z parent population, at least in a statistical 
sense Q. 



1 The assumption that galaxy star-formation and stellar masses 
form a one-parameter family has been cited before in other con- 
texts. It i s implicit, for exampl e, in the star-formi n g main se- 
quence of iNoeske etafl i2007al ) and iNoeske et""aH j2007br i , in 
which galaxies of a given stellar mass today were fitted to a unique 
star-formation hist ory. It is also implic i t in t he halo abundance 
matching model of lConrov fc Wc chslcr (2009), in which they as- 
sumed that the star formation rate was uniquely determined by 
halo mass at each redshift. These and similar prescriptions tend 
to produce a one-parameter family of galaxies labeled by mass. 
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Figure 3. SFR vs. galaxy stellar mass for all model galaxies in the computational box at z = 2 (left), 2 = 3 (centre) and z = 4 (right). 
Each galaxy is shown with a red or a blue symbol depending whether its descendant at z = lies on the red sequence or the blue cloud. 
Symbols surrounded by a black circle correspond to ongoing mergers with mass ratio greater than 1:4 (major mergers). The solid lines 
show, for each redshift, the best linear least-square fit to the main sequence of SFGs, while the dashed lines mark ±3cr. 
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Figure 4. The model results shown in Fig. [3] are c ompared with the galaxies of Santini et al . (2009; z = 2) and different observational 
determinations of the S FR-Af, relati on at z = 2 llDaddi et alj|2007l: lRodighiero et al.ll2010l : iReddv et aTll2012ri , z = 3 l lMagdis et al.l 
l2010l : IWuvts et alJlioilh , and z = 4 jStark et al.ll2009l : iGonzalez et al.ll2012ri . The solid and dotted lines at z = 2 3, 4 show the model 
relation and are the same as in Fig. [3] They have been added to facilitate the comparison with Fig. \3\ 



This kind of diagram is useful to study also the evo- 
lution in the star formation properties of the galaxy popu- 
lation. The tracks in Fig. [2] are the same backward evolu- 
tionary tracks that we already showed in the right panel of 
Fig. [1] but now we have added colour to display the median 
star formation rate (SFR; left panel), the median specific 
star formation rate (SSFR; middle panel), and the fraction 
of red galaxies (right panel) on each mosaic tile. It is impor- 
tant to keep in mind that the tracks on the three panels of 
Fig. [2] are all and only for galaxies that end up on the red 
sequence at 2 = 0. This is why, in the right panel of Fig. [2] 
the fraction of red galaxies at z — is unity for all values of 
stellar mass. 

Here a red galaxy is defined to be in the red sequence 
of the colour-mass diagram. We chose the division between 
the red sequence and blue cloud by eye as the criterion to 
separate red and blue galaxies at each timestep. 

The left panel of Fig. [2]shows that, in our model, there 
is an almost one-to-one relation between the passive popula- 



tion of z ~ massive galaxies (the galaxies with the lowest 
SFRs ever) and the population of massive z = 3 — 6 galax- 
ies with SFR > 100 Mq yr _1 (which comprises the strongest 
starbursts in the Universe). While at low-z, the most mas- 
sive galaxies are those with the lowest SFRs, at z > 3 the 
highest SFRs are predicted to be in the most massive galax- 
ies at those redshifts. 

The right panel of Fig. [5] shows that, in our model, 
the most massive galaxies grow along the blue sequence un- 
til z ~ 3.5, and that the red sequence begins to emerge at 
z ~ 3—3.5. After this point, most giant ellipticals evolve pas- 
sively. One can see clearly from Fig. [2] that the point where 
the SFR drops and galaxy move to the red sequence is ex- 
actly where the curves Af,(z) go horizontal. This is logical, 
since the stellar mass stops increasing when star formation 
shuts down. This shut down is a direct consequence of the 
assumption that gas accretion and star formation shut down 
above a critical mass Merit (z) , which underpin s all t he re- 
sults of ICattaneo et all (|200fj| ), ICattaneo et all (|200Sl ), and 
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this article. Restart of galaxy growth at z ^ 0.5 is due to 
dry mergers at the centres of groups and clusters. This ef- 
fect is likely to be overestimated beca use our model does no t 
include tidal stripping (discussion in lCattaneo et aLll2008l ). 

The middle panel of Fig. [2] shows a trend, in which the 
typical SSFR decreases with time and lower at higher stellar 
masses. 



4 HIGH-REDSHIFT STAR-FORMERS AND 
THEIR DESCENDANTS 

Having described how we map the high-z galaxy popula- 
tion into the low-z one and vice versa, we are now ready 
to present our results. In this Section, we look at high-z 
(z ~ 2 — 4) SFGs and we follow their descendants at z = 0. 

Fig. O shows the SFR - stellar mass relation at z = 
2, 3, 4. Most galaxies (symbols) lie on a diagonal strip where 
the SFR grows with the stellar mass M*. Thi s strip is the 
mam sequence of SFGs (|Noeske et all l2007bl ; lElbaz et all 
l201ll ). Between z = 4 and z = 3, a second strip emerges 
parallel to the former, the tip of which is visible in the panel 
at 2 = 3. It is the sequence of passive galaxies. Galaxies on 
the latter have much lower SFR for a given stellar mass and 
this is where, already at z = 3, the most massive galaxies 
lie. 

We define the SFG population as follows. We sepa- 
rate the main sequence of SFGs from the passive popula- 
tion with the line log SFR = logAf* + c(z) where c(z) = 
(-10, -9.7, -9.5) for z = (2, 3, 4). We then perform a linear 
least squares fit to the galaxies that lie above this divid- 
ing line and that are more massive than 1O 1O M0. The fit 
is the solid line in Fig. [3] We compute the standard devia- 
tion a of the points from the fit and consider the galaxies 
within 3cr of the fit to be the main sequence of SFGs (the 
galaxy population within the two dashed lines). We define 
an SFG to be any galaxy above the lower dashed line, includ- 
ing the outliers from the main sequence of SFGs. We, there- 
fore, do not adopt an absolute SFR criterion to define SFGs, 
though we introduce a stellar-mass cut at M» > 10 10 Mq. 
One can see from Fig. [3] that, with this cut, our simulated- 
galaxy sample is 'complete' down to ~ (100, 60, 30) Mq yr -1 
at z — (4, 3, 2). In the next Section, we shall start from mas- 
sive red galaxies at z — and we shall track their progenitors 
back in time. 

In our model, the normalisation of the SFR - M* rela- 
tion decreases from high to low z for both star-forming and 
passive galaxies. A least squares fit to the main sequence 
of SFGs, LogSFR/(M yr- 1 ) = a{ Log M./M© - 11) + 6 
yields similar slope of a ~ 0.80 in each redshift bin, but 
the normalisation at M» = 10 11 Mq decreases from SFR ~ 
160 M Q yr _1 at 2 = 4 to SFR ~ 100 Mq yr" 1 at z = 3 and 
SFR ~ 40 Mq yr -1 at z = 2 Hence, the number of galaxies 
above a given SFR decreases with time. For instance, the 
number density of galaxies with SFR > 100 Mq yr -1 is pre- 
dicted to decrease from ~ 170 per (150 Mpc) 3 at 2 = 4 to 
~ 50 and ~ 7 per (150 Mpc) 3 at z — 3 and z — 2, respec- 
tively. 

Fig. [4] compares the model results of Fig. [3] with obser- 
vational determinations of the SFR - M* relation. Daddi et 
al. (2007; 2 = 2), Santini et al. (2009; = 2), Rodighiero 
et al. (2010; 2 = 2), Reddy et al. (2012; 2 = 2), Stark et 



al. (2009; z = 4), and Gonzalez et al. (2012; 2 = 4 used 
the Salpeter IMF in their estimates of M*, while Magdis 
et al. (2010; 2 = 3) and Wuyts et al. (2011; z = 3) used 
the Chabrier IMF. We have scaled t heir values to be consis - 
tent with the Kenni c utt IM F used in lCattaneo et all (|2006l ) 
and lCattaneo et all (|2008l ). To pass from the Salpeter IMF 
and the Chabrier IM F to the K e nnicu t t IMF, we have use d 
the relations found in| Bell et al.l (f2003h ; iFardal et all ([2007b : 
iTrever et all (120071 ). 

The slope and the scatter of the model relations at 2 = 2 
and 2 = 3 are quite similar to those that we find in the 
data. However, the normalisation of the SFR - M* relation 
at M, = 1O 11 M0 that we find in our model at 2 = 2 is 
~ 2 — 3 times lower than that found in observational stud- 
ies dDaddi et alj|2007l: ISantini et alj|2009l : iRodighiero et all 
l20ld : iReddv et all 120121 )7 We find the same discrepancy at 
z — 3, when we c ompare our re s ults t o the observational 
determ inations by iMagdis et al.l (|2010T ) and IWuvts et al.l 
|201l| y 

At z = 4 our SFR - M* relation ap pears to be higher 
tha n the observational dete rminations bv lStark et al.l (| 20091 ) 
and iGonzalez et al.l (|2012T l. However, the comparison with 
the data at z = 4 requires extreme caution for two reasons. 
First, the blue and the red hatched ares finish abruptly at 
~ 10 1 ' 5 Mo/yr and ~ 10 1 ' 8 Mq /yr, respectively, because the 
relation at 2 = 4 is given in bins of SFR, rather than in bins 
of M* . Secondly, the data bv lStark et al. I (|2009D are not dust 
corrected. 

Globally, the model is not doing too bad considering 
that the rough agreement in Fig. 0] has been obtained with- 
out tuning any free parameter . We have simply taken a 
model that is already published jCattaneo et alj|2006l . l2008h 
and that works well at z = 0, and we have analysed its pre- 
dictions for the SFR - M* relation at 2 = 2 — 4. Admittedly 
we have also required that this model reproduces the lumi- 
nosity function of Lyman-break galaxies at 1700 angstrom 
restframe. However, this is no guarantee to reproduce the 
SFR - M* relation, particularly since many SFR measure- 
ments come from infrared data. 

In closer detail, the discrepancy with the observed SFR 

- M* relation at 2 = 2 was not unexpected because we know 
that the cosmic SFR density predicted by our model at 2 = 2 
is at the lower limit of the range allowed by observations 
(Fig. 8 of Cattaneo et al. 2006), and that all semianalytic 
models underpredict sub-mm counts, u nless they invoke a 
top-heavy stellar initial mass function (|Baugh et all 120051 : 
lLacev et ai]|2008h . 

We should also note, however , that at 2 = 1, where 
there are more data, iNoeske et al.l |2007bl ) and IChen et al.l 
|2009l ) find a subs tantially lower norm alisation of the SFR 

- M* relation than lElbaz et al.l (|2007l ) and Santini et al. do 
(Elbaz et al. 2007 and Daddi et al. 2007 are the same group). 

Having discussed the comparison with the data, we now 
take a closer look at the results of our model. 

Fig. 3 shows that, at 2 = 4, a lot of the SFGs with SFR> 
100 Mq yr -1 are objects at the massive end of the SFR se- 
quence. In contrast, at 2 = 2, SFGs with SFR> 100 Mq yr" 1 
are outliers from the main sequence. 

The main sequence of SFGs represents the normal 
SFG population at a given redshift. Yet, at 2 ~ 4, the 
tip of the main sequence of SFGs reaches SFRs of about 
1000 M yr -1 . 
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Figure 5. Fraction of galaxies on the red sequence at z = as a 
function of peak star formation rate. The error bars arc Poissonian 
errors. This figure includes all galaxies with M, > 10 10 Mq at 

2 = 0. 



Some galaxies have very high SFRs (up to thousands of 
Solar masses per year) that cause them to be outliers from 
the main sequence of SFGs. It is apparent from Fig. 0that 
these galaxies are almost always mergers, because objects 
which have experienced a greater than 1:4 merger in the 
last 300 Myr have been surrounded with a black circle. 

Fig. suggests that mergers are necessary to explain 
the highest SFRs at any given z, but that they do not drive 
the bulk of star formation at high redshift, most of which 
occurs along the main sequence of SFGs. The circled symbols 
account for 0.6, 1.2, and 1.4% of the objects on the main 
sequence of SFGs (i.e., those between the 3a lines in Fig. [3} 
at z = 2, 3, 4, respectively. The percentages increase to 54, 
92, and 100% if we look at outliers above the main sequence 
of SFGs, defined as galaxies that lie more than 3<r above the 
fit to the main sequence of SFGs (upper line in Fig. [3j ■ 

Galaxies have been plotted in red or blue depending 
on whether their descendants lie on the red sequence or the 
blue cloud of the colour - M» diagram at z = 0. The criterion 
used to separate red and blue galaxies is U — B > 0.95 + 
O.O6(logM J ,/M - 10.0). 

A strong result that emerges from Fig. 3 is that virtually 
all massive SFG galaxies evolve into the red sequence at 
= 0. The fraction of SFGs that end up on the red sequence 
depends slightly on mass, SFR and redshift. We find that 
~ 95 - 98% of high- 2 SFGs with M* > 10 11 Mq evolve into 
red-sequence galaxies by z = (these values correspond to 
z = 2 and z = 4, respectively). At lower masses (10 Mq < 
M, < 10 11 Mq), the fraction of SFGs that end up of the red 
sequence is lower and exhibits a slightly stronger dependence 
on redshift. The SFGs with 10 10 Mq < M» < 10 11 Mq that 
evolve into red-sequence galaxies by z = are 86%, 91%, 
and 97% at z = 2, 3, 4, respectively. 

Above a given SFR, say lOOMQyr -1 , objects have a 
substantial probability of ending up on the red sequence 
at z = (96%, 95% and 97% for objects with SFR> 
100 M Q yr" 1 at z = 2, 3, 4, respectively). 

Fig. 5 shows the probability that a z = galaxy is on 
the red sequence as a function of its peak SFR. It shows 



that 80% of all objects that have ever reached > 40MQyr _1 
evolve into red-sequence galaxies. 

We do not argue that the peak SFR of a galaxy is the 
most important factor in determining if it ends up on the red 
sequence or the blue cloud. But since the SFR of a galaxy at 
any time put a lower limit to its peak SFR, Fig. 4 is useful 
because, given a SFR measurement, it gives the minimum 
probability that an observed galaxy has to end up on the 
red sequence at 2 = 0. 

In fact, we do not argue at all that absolute SFR 
provides a physical criterion to tell if a galaxy becomes 
quenched. Galaxy mass (due to its link to halo mass) and 
morphology are, in this sense, much more relevant. The pur- 
pose of this analysis is rather to show the correlation be- 
tween one of the most basic properties that one can measure 
in a high-2 galaxy (SFR) and what our model predicts for 
its later evolution. 



5 THE PROGENITORS OF MASSIVE RED 
GALAXIES 

In Section 3, we have looked at high-redshift (z ~ 2 — 4) 
SFGs with M„ > 10 10 M and we have found that most 
of them evolve into massive red-sequence galaxies by z = 0. 
We now do the opposite. We start from massive red-sequence 
galaxies in the local Universe and we reconstruct their SFR 
histories by tracking their most massive progenitor back in 
time. 

Fig. [6] has been constructed by selecting eight galaxies 
with M, it 10 11 ' 5 Mq today. These galaxies have been cho- 
sen to illustrate the variety of star formation histories of 
today's massive red galaxies The total SFR as a function 
of 2 (black solid line) has been decomposed into the contri- 
butions of disc and bulge. Mergers with mass ratios greater 
than 1:4 (major mergers) are highlighted with an asterisk 
symbol. The vertical dashed lines show the redshifts at which 
the dark matter haloes of these galaxies have passed M cr i t . 

In galaxies (a), (b) and (c), disc star formation dom- 
inates the total SFR before star formation is quenched at 
3 < 2 < 4 (the fact that the bulge SFR is higher than 
the disc SFR after quenching is an artifact due to assum- 
ing a gas-surface-density threshold for disc star formation 
but not for bulge star formation). In galaxy (d), the two are 
comparable just before star formation is quenched. Galaxy 
(c) experienced no major mergers (in Fig. [6] major mergers 
are marked with asterisks). Galaxies (a), (b) and (d) ex- 
perienced one major merger each at z < 1, long after they 
ceased to make stars. In conclusion, the mergers experienced 
by galaxies (a), (b), (c) and (d) are dissipationless and the 
bulge star formation in galaxies (a), (b), (c) and (d) is driven 
not by mergers but rather by disc instabilities. 

Galaxies (e), (f), (g) and (h) had 2-3 major mergers 
each and they exhibit strong SFR peaks at the time of their 
first merger, which occurs before quenching, when the galax- 
ies were still gas-rich. There is a clear direct link between 
dissipative mergers and strong starbursts in the bulge com- 
ponent, which cause it to dominate the SFR temporarily. In 
galaxies (e) and (f ) , mergers cause the SFR to increase by a 
factor of ~ 4. In galaxy (h), the SFR increases by a factor of 
~ 10. It is also interesting to note that while galaxies (e), (f) 
and (g) go back to normal 'quiescent' star formation after a 
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Figure 6. The SFR histories of eight red-sequence galaxies selected for illustrative purposes. Each SFR history has been decomposed 
into disc (blue) and bulge (red) star formation. The stellar mass of each galaxy at 2 = has been shown in the corresponding panel. 
Major mergers (with ratio > 1 : 4) have been marked with asterisks. Most of them are 'dry' (dissipationlcss) mergers occurring after the 
galaxies have ceased to make stars. 'Wet' (gas-rich) mergers are accompanied by strong starbursts, characterised by a sharp rise and a 
fall in the bulge SFR. The fact that the bulge SFR declines continuously and is higher than the disc SFR after quenching is an artifact 
due to assuming a gas-surfacc-dcnsity threshold for disc star formation but not for bulge star formation, which, after shutdown, is fed 
by stellar mass loss only. 



merger-driven starburst - if one can call quiescent a galaxy 
like (g) with a disc SFR > 100 M© yr" 1 ), - in galaxy (h) the 
merger-driven starburst coincides with the shutdown of star 
formation. 

Merger-driven starbursts are, by construction, short- 
lived. We have verified that in galaxies (e), (f) and (g) only 
a small fraction of the final galaxy stellar mass was formed 
in the bulge even if the highest SFR occurred in the bulge. 
Galaxy (h) is the prototypical example of an elliptical galaxy 
that makes most of its stars in a single merger-driven burst 
at high z. However, objects of this type are exceedingly rare. 
They can be counted on the fingers of one hand in our com- 
putational volume of (150 Mpc) 3 . 

With insight from these particular cases, we are now 
ready to look at the statistical properties of the progeni- 
tors of red-sequence galaxies. We split the latter into three 
ranges of stellar mass at z = (lower intermediate mass: 
1O 1O ' 5 M < M, < IQ^Mq; upper intermediate mass: 
1O U M < M, < 1O 11S M ; high mass: h'L > 1O U - 5 M ) 
and characterize the SFR histories of the disc and the bulge 
of each galaxy in terms of two properties: the peak SFR and 
the redshift at which the SFR reaches its peak value. 

Fig.Qshows the joint distribution of SFR pca k and z pea k 
for galaxies in which in the peak SFR is linked to star forma- 
tion in the bulge component (above) and to the disc compo- 
nent (below). The histograms in Figs. [8][9] show the distribu- 
tions for SFRp Ca k and z pea k separately. In these figures, blue 
histograms are those galaxies whose peak SFR was in the 
disc (they correspond to the lower panels of Fig. [7]), while 



red histograms are those whose peak SFR was in the bulge 
(they correspond to the upper panels of Fig. [7]). 

Figs. 17191 show that the typical peak SFR is tens of So- 
lar masses per year for lower intermediate-mass red galax- 
ies, about a hundred Solar masses per year for upper 
intermediate-mass red galaxies and several hundreds Solar 
masses per year for high-mass galaxies, but this is only a 
statistical statement. Even among lower intermediate-mass 
galaxies, there are objects with SFR pca k ~ 1000 M yr -1 . In 
all ranges of mass, galaxies in which the peak SFR occurs in 
the disc are much more numerous than galaxies in which the 
peak SFR occurs in the bulge. To put it quantitatively, only 
11%, 12%, and 16% of the red-sequence galaxies in the stel- 
lar mass intervals 10 10 ' 5 M < M» < 10 11 M Q , 10 11 M < 
M» < 10 11 ' 5 M Q , and 10 11 ' 5 M Q < M* < 10 12 M have ex- 
perienced a SFR peak dominated by bulge star formation. 

As we argued at the end of Section 2.3, disc star for- 
mation is always stream-fed in the sense that, by construc- 
tion, discs can acquire gas from smooth accretion only. In 
contrast, both mergers and rapid smooth gas accretion con- 
tribute to bulge star formation. Major mergers do it directly, 
by bringing gas into the starburst /bulge component. Gas ac- 
cretion and minor mergers do it indirectly, by causing a sud- 
den increase in the gas content of the disc, which becomes 
unstable and transfers mass to the bulge. Thus the counts in 
the upper panel of Fig. [7]and the red histograms in Figs. 18191 
represent upper limits to merger-driven SFR, and we refer to 
these histograms loosely as the merger-driven contribution. 

Therefore, our model makes two strong predictions. The 
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Figure 7. The joint distribution of peak SFR and rcdshift of peak SFR for z = red-sequence galaxies in three bins of z = stellar 
mass: 10 10 - 5 M Q < M* < lO n M , 10 11 Mq < Af* < 10 115 M Q , and M* > lO u - B Af . The upper and lower panels show the SFR pcak - 
z pcak distribution for galaxies with SFR^fJ^ > SFR^J eafc and SFR^f^ < SFR^ , respectively. 
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Figure 8. Distribution for the value of the peak SFR for red-sequence galaxies in three intervals of stellar mass at z = 0: 10 10 5 Mq < 
M* < 10 11 Mq (left), 10 11 M Q < M» < 10 11 ' 5 Mq (centre) and 10 11 - 5 Mq < M, < 10 12 Mq (right). The distribution of SFFi pcak has 
been shown separately for disc- dominated peaks (blue histograms) and bulge-dominated peaks (red histograms). 



Galaxies on the Red Sequence Today 




Figure 9. Distribution for the value of the redshift Zpeak at which the SFR reaches its maximum for red-sequence galaxies in three intervals 
of stellar mass at z = 0: 10 10 ' 5 Mq < M» < 10 11 Mq (left), 10 11 Mq < M» < 10 11 - 5 Mq (centre) and 10 11 ' 5 Mq < M« < 10 12 Mq 
(right). The distribution of 2 pca ]j has been shown separately for disc-dominated peaks (blue histograms) and bulge-dominated peaks (red 
histograms). 
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Figure 10. The contribution of bulge star formation to the total in situ star formation of red sequence galaxies (computed by following 
the galaxies' main progenitors) . The distribution of values for J SFRi) U ig e (t) At/ J SFRtot (t) At is shown for three bins of stellar mass 
at z = 0: 10 10 - 5 Af© < Ah < 10 11 M Q , 10 11 M© < Af» < 10 11 - 5 M e , and M, > 10 11 - 5 Af©. 



first is that, in most red galaxies, star formation was stream- 
fed, not merger-driven. The second is that, conversely, merg- 
ers are responsible for the most intense episodes of star for- 
mation at each redshift see Fig. [8] in agreement with data 
over a broad range of redshifts from z ~ to z ~ 2. 

The first conclusion can also be verified directly by tak- 
ing the total SFR history and the bulge SFR history of each 
galaxy, computed by tracking the main progenitor back in 
time (i.e. the equivalent of the black and the red curves in 
Fig. [6]), and by integrating them over time, to determine 
what fraction of the star formation that has occurred has 
been bulge star formation. Fig. 1 101 shows the distribution of 
values for f*° SFR bu i gc (t) At/ f*° SFR tot (t)d£ in the three 
bins of galaxy stellar mass 10 10 ' 5 M© < M* < TO 11 M©, 
10 n M© < M, < 10 n ' 5 M©, M, > 10 n ' 5 M© (stellar 
masses at z = 0; to is the current age of the Universe; 
and this figure is only for red-sequence galaxies). The av- 
erage contribution of merger driven star formation for the 
three stellar mass bins is 11%, 15%, and 26%, respectively. 
Galaxies that formed most of their stars in merger-driven 
starbursts are rare, though it is quite usual that a quarter 
of the in-situ star formation in a massive early-type galaxy 
occurs via this mode. Most of the stars in a giant elliptical 
were not formed in situ, i.e. they were formed in smaller 
objects that merged with the galaxy, so we do not include 
either in the ~ 1/4 of bulge star formation or the ~ 3/4 of 
disc star formation. 

Fig. [5] shows the distribution for the redshift z pca k at 
which the peak SFR occurs for the three ranges of masses. 
For lower intermediate-mass red galaxies, z pca k ~ 2. For 



upper intermediate-mass red galaxies, z pca k ~ 2.5 — 3.5. For 
high-mass red galaxies, £ pea k ~ 4. 

The decrease of z pea k at lower masses is an aspect of 
downsizing. In our model, downsizing occurs because haloes 
with a lower mass at fixed z cross M cr i t at a lower redshift 
(Fig. 4 of Cattaneo et al. 2008). In the same way M cr i t is 
constant at z 3 and increases at higher z, this also applies 
to the characteristic stellar mass M" 1 with which galaxies 
enter the red sequence. The characteristic star formation 
rate of a galaxy that enters the red sequence at cosmic time 
t is SFR ~ Merit (t)/t, which must decrease with t, since 
~ M cr it (t) is a non-growing function of t. 

In lCattaneo et ail (|2008T ). we analysed the time of peak 
SFR as a function of the final stellar mass for massive red- 
sequence galaxies. The results w ere found to be in agreement 
with that iThomas et all (|2005l ) inferred from the spectra of 
local early- type galaxies (Figs. 7-8 of Cattaneo et al. 2008). 

Stream-fed star formation is responsible for most of the 
SFR peaks at all redshifts except perhaps at low-redshift 
(z < 1). In the stellar mass bins 10 10 S M© < M, < 
lO^M©, lO^M© < M„ < 10 n - 5 M Q , and 10 11 ' 5 M© < 
M* < 10 12 M©, the red galaxies with 2 pcak < 1 are 9%, 5%, 
and 2%, respectively. If we look at this specific population, 
then the percentages of galaxies with merger-driven SFR 
peaks rise to 56%, 80% and 67% (instead of 11%, 12%, and 
16%, as we saw before). These values would have been even 
higher if we took z < 0.5 instead of z < 1 (we would have 
found 81%, 98% and 100%, respectively). Fig.[7]tells us that 
these merger-driven peaks are typically about 50M©yr _1 
for the lower-intermediate mass bin, and about 100M©yr _1 
for the upper intermediate mass bin. Merger-driven star 
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formation is therefore relevant to those M* ~ Mq el- 
lipticals that have formed via low-redshift gas-rich merg- 
ers such as those th at we witness in local ULIRGs (see 
ICattaneo et al.ll201ll for a discussion of the role of gas-rich 
mergers in the build-up of the galaxy population) . The same 
trend is also seen in massive red-sequence galaxies: the SFR 
peak is usually stream-fed and is usually at z ~ 4 but in 
the few objects with merger-driven SFR peaks, the SFR 
peak is normally around z ~ 3. Most of the merger activ- 
ity of massive red galaxies is predicte d to be dissipationless 
and to occur at fairly lo w redshifts (|Cattaneo et alJl200St 
iDe Lucia fc Blaizodl2007h . 



6 DISCUSSION AND CONCLUSION 

To summarize, we find that: i) most galaxies that have 
experienced SFRs > 40-MQyr -1 evolve into red-sequence 
galaxies, ii) 80% of the red-sequence galaxies with M* > 
10 11 M Q passed through a phase of high SFR (SFR peak > 
40 Mq yr _1 ) at z > 1 , and iii) in ^ 90% of these cases, this 
phase was stream-fed (merger ratio at peak SFR less than 
1:4). 

Still, mergers are necessary to explain the galaxies with 
the highest SFRs (e.g., the rare objects with SFRs of thou- 
sands of Solar masses per year). Furthermore, if we con- 
centrate on z — red-sequence galaxies with peak SFR at 
2 < 1 (9%, 5%, and 2% of the objects in the three mass in- 
tervals), then the percentages of merger-driven SFR peaks 
rise to 56%, 80% and 67%, respectively. These values would 
have been even higher if we took z < 0.5 instead of z < 1. 
Hence, mergers are important for that small percentage of 
latter-day ellipticals that form in local ULIRGs, even though 
this cannot be the path via which most ellipticals formed 
jOstrikdll98Cl ). 

These figures (and all the others in this article) should 
be taken with caution because our model contains many 
arbitrary assumptions. Their uncertainties cannot be eas- 
ily quantified but we believe they are probably greater than 
10%. The purpose of the percentages quoted in this article 
is to illustrate the magnitude of the discussed phenomena. 
A change from 56% to 67% in the fraction of merger-driven 
peaks from one mass bin to another is not significant when 
one considers the uncertainties of the model. However, we 
prefer to give the exact numbers returned by our model ac- 
companied by this disclaimer than to use vague language, 
where any quantitative information is lost. 

The accuracy of a semianalytic model is established a 
posteriori by how well it the reproduces the observations. 
Our model is in very good agreement with local data. The 
main problem is the discrepancy with the normalisation of 
the SFR - M* relation at z = 2 — 3, which we have discussed 
in Section 4. 

At M* = 10 u M o , our model predicts SFRs that 
are ~ 2 — 3 times lowe r than those inferred obser - 
vationally by iDaddi et all J2007t). ISantini et all J2009I) 



found in the Munich model (e.g. iKitzbichler fc WhitdkoOTT l 
and in the Durham model, which proposes a top-heavy stel- 
lar initial m ass function 
iBaugh et al.l 



120051 : [Lacev 



., . . — ~~~~ . J ~~~~ 

(IMF) as a possible solution (e.g. 
et al.ll200Sl). 



iRodighiero et all (|2010D. iMagdis et all (|2010l ). IWuvts et all 
( 201 lft . lReddv et all (|20ia v" and, at the massive end of the 
mam sequence of SFGs, by iHuang etafl (|2009r i. Our av- 
erage SFR at a given stellar mass would still be too low 
compared to these data even if all the accreted gas were im- 
mediately converted into stars. The same basic problem is 



Since iNoeskeet all (|2007al ) and I Chen etafl (|2009l ) find 
a lower normalisation at z — 1 than the authors mentioned 
above, it is possible that these authors' SFR measurements 
at 2 = 2 need to be lowered, too, or that their masses 
are systematically underestimated. However, we also know 
that semianalytic models underpredict sub-mm counts, un- 
less they invoke a top-heavy IMF. So either a top-heavy 
IMF is the solution, in which case there may be nothing 
wrong with our predicted SFRs, and it is the observation- 
ally derived SFRs that need to be corrected accordingly, 
or this signals a problem in our galaxy formation model. 
For instance, it is possible that our specific SFRs at z = 4 
are too high (Fig. 4) because gas accre ted at z ^ 4 is con- 
verte d into stars too rapidly (see, e.g., iKrumholz fc Dekell 
2012). and this may be the reason why, in the model, there 
is not enough star formation at later epochs (z ~ 2 — 3; 
Fig. [3) , though it also possible that 2 = 4 SFRs inferred by 
IStarket all (|2009l ') are underestimated due to dust. 

In fact, one should also remember that the observa- 
tional determinations of the SFR-M* relation are themselves 
model-dependent. Measured SFRs and stellar masses de- 
pend on stellar-population synthesis models, dust models, 
and assumptions about the galaxies' SFR histories. 

The question here is: to what extent these open issues 
affect our confidence in our results? 

Could the underpredicted SFRs for a given M* be due 
to an underestimate of the importance of mergers? Could 
our conclusion that the bulk of star formation in the ances- 
tors of red galaxies is stream-fed be a result of our model 
somehow overestimating the importance of the stream-fed 
mode? 

It is not inconceivable that our model may somehow un- 
derestimate the merger rate at high redshift, even though the 
assumptions that it makes are quite standard. The problem 
is that the entire main sequence of SFGs is shifted with re- 
spect to the data at z = 2. Mergers can temporarily increase 
the SFR by accelerating the conversion of gas into stars but 
they cannot change the average SFR for a given galaxy mass 
because a shorter star formation timescale means that at the 
end less fuel for star formation is left. 

Let us now consider the second possibility, i.e. that the 
importance of stream-fed star formation is overestimated. 
Even though our model incorporates disc instabilities, in 
this article we have effectively identified stream-fed star 
formation with disc star formation and merger-driven star 
formation with bulge star formation. In reality, while disc 
star formation is always stream-fed, bulge star formation 
is partly merger-driven and partly stream-fed. By equat- 
ing stream-fed star formation with disc star formation and 
merger-driven star formation with bulge star formation, we 
should underestimate the former and overestimate the lat- 
ter, not the opposite, but let us analyse this point in closer 
detail. 

Our semianalytic model has been run using the disc 
stability criterion Q for a disc instability threshold param- 
eter of r\ = 0.7. In other words, we assume that, within 
the disc half-mass radius, the stea dy state is that the mass 
in the disc is half the total mass. iDekel et al.1 (2009) have 
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Figure 11. The gas content of major mergers at z = 2 3, 4. GallCS galaxies that have experienced a merger in the last 300 Myr are 
shown as circles on the SFR- A'/* diagram. Mergers observed in Tacconi et aL, (2008) are marked as stars (six-point stars correspond to 
upper limits on the gas fraction). Their stellar masses and SFRs have been adjusted to the Kennicutt IMF. Each symbol is colour-coded 
according to the merger gas content, i.e., M gas /(M, + M gas ). jwbA circled star symbol indicates that the gas fraction is an upper limit, 
while the rest of the data from Tacconi et al. have uncertainties of about 0.08 in the fraction. The lines show the position of the main 
sequence of SF galaxies. They are the same as in Figs. \3\ and [4] 



explored a model for the instability of clumpy discs and 
they have been able to relate the disc-to-total mass ratio 5 
within the the disc half-mass radius to a/v Iot , the ratio of the 
velocity dispersion a within the disc to the disc's rotation 
speed w r ot (also see ICenel et alj|2012h . For a Toomre disc 



1964 



Tremainel 

1 and a flat rotation curve, their 



Binncv 



insta bility parameter Q (|Toomri 
|200S| , Chapter 6) of Q 
model gives 8 ~ \/2 ■ a/v TOt . So, rj = 0.7 is equivalent to 
cr/urot — 0.35. In fact, spectroscopic observations of gas 
kinematics in massive discs at z ~ 2 give values closer to 
a/vrot ^ 0.2 (lErb et al.ll2004l ; iForster Schreiber et al.ll200d : 
ICresci et al.ll2009h . which imply S ~ 0.3. Based on this ana- 
lyis, the actual inflow due to disk instability in high-z clumpy 
discs is stronger than what GallCS assumes. 

Still, our predicted bulge-to-total stellar mass ratios 
for blue galaxies are systematically on the high side at all 
masses (GallCS predicts a typical bulge-to-total mass ra- 
tio of Mbuige/Afgai ~ 0.3 fo r a blue galaxy with M r ~ 
-20.5; ICattaneo etall 120061 . Fig. 10). Simply lowering 
the value of 77 would bring predicted bulge-to-total stel- 
lar mass ratios to values that are incompatible with the 
prevalence of late-type morphologies in the local galaxy 
population. This is why hydrodynamic simulations have 
for a long time experienced major difficulty in forming 
spirals with acceptable bulge-to-total stellar mass ratios, 
though they follow gas and stellar dynamics more self- 
consistently than semianalytic models do. Recent progress 
in forming spirals with acceptable bulge-to-total stellar 
mass ratios is linked to strong feedb ack that preferentially 
ejects gas from the central starburst (Governato et all 20ld; 



Guedes et al. 201ll;|Piontek fc Steinmetal201ll ; lBrown et al.l 
20071 ; iMcCarthv et al.ll2012l ). Therefore, a more physical de- 



scription of disc instability must be accompanied by a more 
physical description of stellar feedback. 

Based on this discussion, GallCS likely underestimates 
both the inflow into the bulge due to disc instability and 
the outflow from the bulge due to stellar feedback. So, part 
of the stream-fed star formation that we predict to occur in 
the disc may have occurred in the bulge, or vice versa (as 



it is more likely, at least for spiral galaxies, whose bulge-to- 
stellar mass ratio are overpredicted). However, that does not 
change the relative importance of stream-fed and merger- 
driven star formation. 

One may also worry that our SFRs have been boosted 
by a factor of (1 + z) ' 6 to produce m ore star formation at 
high redshift (Eq. 2; also see Fig. 8 of lCattaneo et al]|2006l 
for the impact of this assumption on the evolution of the 
cosmic SFR density). While this assumption increases the 
stream-fed SFR, the factor (1 + z) ' 6 has been assumed to 
multiply both the disc SFR and the bulge SFR. Thus it does 
not change the relative importance of the two in terms of 
star formation efficiency. 

The most serious concern is that our model may exag- 
gerate gas consumption in protogalaxies at z ^ 4. with the 
consequence that when the first major mergers occur, there 
is little gas left to trigger a starburst. To address this point, 
we have analysed the gas content of major mergers. 

Fig. [11] is a remake of Fig. [3] where we only show the 
galaxies that have experienced a major merger in the last 
300 Myr (the circled ones) and compare them to the merger- 
like objects observed in Tacconi et al., (2008). We colour- 
code the symbols according to their gas fraction. In GallCS, 
mergers of quenched, red-sequence galaxies are gas poor, but 
that is due to our quenching criterion, which is related to 
the critical halo mass M cr it, and has nothing to with the star 
formation efficiency. So, let us focus on the galaxies that are 
on the main sequence of star forming galaxies. 

In our model, at z = 2 and M* ,$ 10 10 ' 5 M©, gas frac- 
tions range from ~ 0.3 to ~ 0.6 There is a declining trend 
with mass, which is also observed (Tacconi et al. 2012). At 
M» ~ 10 11 M Q , the typical gas fraction is ~ 0.15 - 0.2, in 
agreement wit h the gas fractions of major mergers observed 
at z — 2 — 3 (|Tacconi et al.l [20081 ). Therefore, there is no 
reason to believe that we underestimate the importance of 
merger-driven star formation because the gas fractions of 
our mergers are systematically lower than they should be. 

Our resu lts axe consistent with those that 

ICattaneo et al.l (|201ll ) obtained with a much simpler 
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galaxy formation model. In that article we concentrated 
on mass assembly (here we concentrate on star forma- 
tion). We found that gas-rich mergers make a negligible 
contribution to the baryonic mass assembly of the overall 
galaxy population, although they contribute about half 
of the mass in ellipticals with M* < lO 10,8 ^ -1 M Q , which 
formed their stars at lower z th an giant ellipticals due to 
downsizing (|Thomas et al.l 120051 ') . Elliptical galaxies with 
Alt > lO 11 ft _1 M0 accrete most of their mass via dry 
mergers. In Cattaneo et al. (2011), we did not consider the 
possibility that the critical mass above which gas accretion 
is shut down may be higher at higher redshifts and we 
wondered to what extent that could affect our conclusions. 
Here we have included a redshift dependence of M cr it 
(Eq. 1), but the basic picture has not changed. 

Lowering z cr i t , which is equivalent to increasing M cr it at 
high redshift, delays the shutdown of gas accretion but also 
transforms a number of dry mergers into wet mergers. The 
total amount of star formation increases but the relative im- 
portance of the two modes is almost unchanged. Values of 
Zcrit substantially lower than z ~ 3 cause GallCS to over- 
predict the galaxy luminosity function at z — 0. This effect 
can be partially compensated by lowering the value of M cr it 
at z — but there is limited leeway to do so. 

In conclusion, the parameters of our model contain 
great uncertainties (particularly those anchored to the lumi- 
nosity function of Lyman-break galaxies), which affect our 
capacity to trust the model's predictions in quantitative de- 
tail, especially since we know that it fails to reproduce the 
SFR-Ai, relation inferred from observations at z ~ 1 — 4. 

However, the basic results of our model are quite robust 
to changes in parameter values and model assumptions be- 
cause it is difficult to increase the importance of wet mergers 
without also increasing the importance of stream-fed star 
formation, unless one finds a way to accrete large masses 
of gas while preventing them from making stars until the 
first mergers occur. These result s corroborate the scenario 
portrayed by iDekel et al.l (|2009h . where they made a dis- 
tinction between the overall SFG population and extreme 
SMGs, which are about ten times less numerous. They con- 
clude that the accretion of cold gas is the main mode of 
galaxy formation in the former population, whereas mergers 
are only necessary to explain the latter. 
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